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Abstract: 

Jet event rates in deep inelastic ep scattering at HERA are investigated applying 
the modified JADE jet algorithm. The data are corrected for detector and hadroniza- 
tion effects and then compared with perturbative QCD predictions using next-to- 
leading order calculations. The strong coupling constant a s (M|) is determined eval- 
uating the jet event rates. Values of a s (Q 2 ) are extracted in four different bins of the 
negative squared momentum transfer Q 2 in the range from 40 GeV 2 to 4000 GeV 2 . 
A combined fit of the renormalization group equation to these several a s (Q 2 ) values 
results in a s {M 2 z ) = 0.117 ± 0.003 (stat) ± gloSa ( S V S ) + °- 006 U et algorithm). 
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1 Introduction 



Multi-jet production in neutral current deep inelastic positron proton scattering (DIS) is in- 
vestigated using data taken in 1994 and 1995 with the HI detector at HERA at DESY, where 
820 GeV protons collide with 27.5 GeV positrons. A measurement of the relative production 
rates of events with different jet multiplicities will lead to a better understanding of the underly- 
ing elementary processes sensitive to perturbative QCD (pQCD) including its basic parameter, 
the strong coupling constant. 

The DIS reaction can be written as e+p — * e+X, where X denotes an arbitrary hadronic final 
state. The process is sketched in figure |], in which the four-momentum vectors of the particles 
are given in brackets to define the kinematic quantities used later in the text. The incident 
positron is scattered over a wide range of the negative squared four-momentum transfer Q 2 = 
—q 2 = — (k — k') 2 . The absorption of the virtual photon by the proton produces a hadronic final 
state with a squared invariant mass of W 2 = (q + P) 2 . 




Figure 1: Sketch of neutral current (NC) deep inelastic ep scattering (DIS) reaction ep — > eX, 
where the four-momentum vectors are given in brackets. 



In the range of Q 2 investigated in this paper, neutral current DIS reactions can be described 
as purely electromagnetic scattering of the positron off a quark via single photon exchange^. 
To lowest order, this leads to (1+1) jet events. One jet is caused by the scattered quark, the 
other (+1) by the remaining part of the proton, the proton remnant, which can only be partially 
observed in the HI detector. 

Higher jet multiplicities are described in the framework of perturbative QCD. The corre- 
sponding cross sections are expressed as power series in the strong coupling constant a s , such 
that the leading order (LO) contribution to (2+1) jet events is of order 0(a s ). Two generic 
graphs contribute (figure ||), where the hard scattering process is either gluon (2a) or quark (2b) 
induced. The resulting jets from the hard scattering are called "current jets" . 

By applying the modified JADE jet algorithm (jj ^] to the data it has been shown in previous 
publications || f|, that jet production at HERA is dominated by the two diagrams of fig. || in 
a large region of phase space. This motivates a quantitative comparison to pQCD predictions 
including real and virtual corrections to the diagrams in question. 

Up to now the theoretical calculations || |6| available for such comparisons neglected terms 
containing the product of W 2 and the JADE jet resolution parameter y cu t and thus were not 
valid for large values of W 2 accessed at low values of Q 2 . In recent NLO QCD-calculations @, ||] 
no such terms are neglected. Hence it is appropriate to compare more measurements with these 
new predictions in a larger region of phase space, especially for low Q 2 where higher statistics 
are available. Thus this analysis supersedes the previous HI publication (3|. 

1 In the Q 2 range under study (Q 2 < 4000 GeV 2 with the present statistics) effects due to the exchange of the 
intermediate vector boson Z° can be neglected. 
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Figure 2: Generic diagrams of the QCD LO contributions to (2+1) jet events in neutral current 
DIS. 



A necessary prerequisite to the analysis reported here is to define a phase space region in 
which NLO calculations are able to describe the data after correcting for hadronization and 
detector effects; the dependence of the (2+1) jet event rate on Q 2 , together with other hadronic 
control variables, should be well predicted by the fixed order calculation. 

The presented quantitative pQCD analysis is based on the previous HI publication Q. The 
measured (2+1) jet event rates in different bins of Q 2 are corrected for detector and hadronization 
effects to the parton level. Taking the known parton densities in the proton, values of a s are 
extracted as a function of Q 2 which fit best the corrected jet event rates. It is thus possible 
to investigate the scale dependence of the strong coupling constant using only one observable, 
the (2+1) jet event rate, in a single experiment. Then a s (M|) is determined by performing a 
fit of the QCD parameter A in the renormalization group equation to the independent a s (Q 2 ) 
values. Here Mz is the mass of the neutral vector boson Z°. The use of a recent NLO QCD 
calculation program allows a study of the dependence of the obtained a s (M|) on the choice of 
the jet algorithm. 



2 HI Apparatus 

A detailed discussion of the HI apparatus can be found elsewhere Q. Here emphasis is put on 
describing the main features of those detector relevant to this analysis, which makes use mainly 
of the calorimeters and, to a lesser extent, of the central and backward tracking systems. The 
polar angle, 8, is defined with respect to the proton beam direction, which is identified with the 
+z-axis. 

The inner part of the detector consists of the central tracking chamber supplemented by 
a forward tracking detector and a backward proportional chamber (BPC), covering the polar 
angular ranges 25° < 9 < 155°, 7° < 6 < 25° and 155° < 6 < 175°, respectively. The central 
and forward tracking devices are used to determine the vertex position. The central tracking 
chamber and the BPC, together with the vertex, are used to measure the positron scattering 
angle in the central and backward region. The angular resolution achieved is better than 1 mrad. 

The energy of the scattered positron is measured by the liquid argon (LAr) calorimeter and 
in a backward electromagnetic lead-scintillator calorimeter (BEMC); the LAr calorimeter is also 
used for a measurement of the hadronic energy flow. A superconducting solenoid outside the 
LAr calorimeter provides a uniform magnetic field of 1.15 T parallel to the proton beam axis 
throughout the tracking region. 
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The LAr calorimeter [10] extends over the polar angular range 4° < 9 < 154° with full 
azimuthal coverage. The scattered positron enters the LAr for events in which Q 2 is larger than 
about 100 GeV 2 . The calorimeter consists of an electromagnetic section with lead absorbers, 
corresponding to a depth between 20 and 30 radiation lengths, and a hadronic section with steel 
absorbers. The total depth of the LAr calorimeter varies between 4.5 and 8 hadronic interaction 
lengths. Test beam measurements of LAr calorimeter modules have demonstrated energy resolu- 
tions of a{E)/E « 0.12/ y/Ej GeVeO.Ol for positrons Q and a{E)/E « 0.5/y/E/ GeV0O.O2 
for charged pions [^]. The uncertainty of the absolute energy scale for positrons is 3% deter- 



mined with studies based on the double angle method of kinematic reconstruction [13]. The 
hadronic energy scale has been verified from the balance of transverse momentum between the 
hadronic final state and the scattered positron in DIS events and is known to a precision of 

4% m 



The BEMC |T5[, with a thickness of 21.7 radiation lengths, covers the backward region of 
the detector, 151° < 6 < 176°. It is mainly used to measure positrons from DIS processes at low 
Q 2 . The acceptance region corresponds to Q 2 values in the approximate range 5 < Q 2 < 100 
GeV 2 . A resolution of a(E)/E us 0.10/ y/E/ GeV 0.02 has been achieved. By adjusting the 
measured positron energy spectrum to the kinematic peak the BEMC energy scale is known 
to an accuracy of 1% []l5|| . Since 1995 the BEMC and BPC in the backward region have been 
replaced by a lead/scintillating fibre calorimeter (SPACAL) [jj~6|| and a drift chamber (BDC) pTj ], 



3 Event Selection 

The data are divided into two sub-samples depending on the location of the detected positron 
in the BEMC or in the LAr calorimeter. The scattered positron is defined as either the most 
energetic cluster detected in the BEMC or a contiguous cluster in the LAr whose energy in 
the electromagnetic section exceeds 80% of the total cluster energy. The used data samples 
correspond to integrated luminosities of 2.8 pb _1 taken in 1994 for Q 2 < 100 GeV 2 (positrons in 
BEMC) and of 7.1 pb" 1 taken in 1994 and 1995 for Q 2 > 100 GeV 2 (positron in LAr calorimeter). 
All further selection criteria follow closely those used in the HI measurement of the proton 



structure function F2 [18|. 



Events in the first sub-sample (BEMC) must satisfy the following requirements: 

1. The four- momentum transfer Q 2 must be between 10 and 100 GeV 2 . 

2. The energy E e of the scattered positron must be greater than 14 GeV, corresponding to 
a fractional positron energy loss in the proton rest-frame of y = qP/kP £0.5 (see fig. |l] 
for the meaning of four-momentum vectors), thus eliminating background from photo- 
production, where a hadron is misidentified as a positron, to a negligible level (<1%). In 
addition this cut suppresses radiative DIS events. 

3. The angle of the scattered positron must be in the range 160° < 8 e < 173° to ensure that 
it is well contained in the BEMC. 

Events in the second sub-sample (LAr) must satisfy the following requirements: 

1. Q 2 must be greater than 100 GeV 2 . 

2. The energy of the scattered positron E' e must be greater than 11 GeV, y less than 0.7, 
and the quantity 5 = J2dusters(E — P z ) (with E and P z being the energy and longitudinal 
momentum of the calorimetric energy clusters) is required to be within 38 GeV < 5 < 70 
GeV. These cuts exclude efficiently photo-production events (<1%), in which the scattered 
positron is lost down the beam pipe. As for the cut in y discussed above, the lower <5-cut 
further reduces radiative DIS events. 
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3. The angle 6 e should fulfill 10° < 9 e < 150° to ensure that the positron is contained in the 
LAr calorimeter and to avoid the transition region between the LAr and the BEMC. 

Two additional requirements were imposed on both sub-samples: 

1. An event vertex from charged tracks within ±30 cm around the maximum, z mean , of the 
vertex distribution is required; z mean was located at z = +5 cm (+1 cm) in 1994 (1995). 

2. W 2 > 5000 GeV 2 , where W 2 is calculated using the double angle method. This method 
relies on the angle of the scattered positron and an inclusive angle of the hadronic final 
state, and is independent of the jet classification. The W 2 cut ensures enough hadronic 
activity in the detector and is equivalent to a cut on the lepton variable y (y > 0.055). 

In both samples Q 2 and y are calculated from the energy and angle of the scattered positron. 
The selected event sample contains 112985 events with Q 2 < 100 GeV 2 and 5410 (8280) events 
with 100 GeV 2 < Q 2 < 4000 GeV 2 selected from the data taken in 1994 (1995). 



4 Jet Reconstruction 



All events are subject to a jet classification using the energy clusters in the LAr calorimeter 
always excluding the scattered positron. In order to compare the measurement to QCD pre- 
dictions the same jet classification is applied to Monte-Carlo events after passing the full HI 
detector simulation and after using the same reconstruction algorithm as for the data. The 
influence of hadronization and detector effects on the number of reconstructed jets at this "de- 
tector level" is studied by determining in addition the jets at the parton and the hadron level. 
At the "parton level" the coloured partons generated according to pQCD calculations will form 
the jets, whereas at the "hadron level", the fully generated hadronic final state is analysed. 

The analysis presented here is based on a modified JADE jet algorithm Q. In order to take 
into account the invisible part of the proton remnant r the modified JADE jet algorithm includes 
a pseudo-particle in the clustering procedure representing the missing longitudinal momentum 
in the event 0. The JADE algorithm is applied in the laboratory frame. Two objects i,j with 
a squared mass mfj are combined to one new object if = m 2 j/M 2 e j < y cu t, where y cu t is 
the resolution parameter. The scale M re f is taken to be the invariant mass W of the hadronic 
system X, where W is calculated as the invariant mass of the four- vector sum of all hadronic 
objects. From studies of the resolution parameter [19|, 20], y cut is set to 0.02. This choice ensures 



satisfactory correlations between jets reconstructed at the parton and the detector level, and 
high statistics of (2+1) jet events, keeping the contribution of (3+1) jet events below 4%. 

Besides the modified JADE jet algorithm two of its variants, the E0 and the P algorithm, are 
investigated, which differ in the combination of two objects i and j to a new object k. For the 
JADE algorithm the new four-momentum vector p^ is given by Pk = Pi+ Pj , whereas for the E0 
and P algorithm the new four-momenta are massless and given by p^ = [Ei+Ej, (Ei+Ej) ■ i ] 
and pk = [\pi +Pj\,Pi + Pj]i respectively. The resolution criterion for all of these algorithms is 
given by yij = 2EiEj(l — cos 9ij)/W 2 , where E^,Ej denote the energies and d%j is the angle 
between the two objects i,j under consideration. As long as there is no restriction on the 
momenta of the resulting jets, the JADE and E0 algorithm are equivalent in the theoretical 
NLO (2+1) jet cross section calculations. For these two algorithms, only the (1+1) jet cross 
section is affected differently by cuts on the current jet angles. The difference, which is marginal, 
is caused by events where the clustering of all three partons and the proton remnant leads to a 
(1+1) jet event. 
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5 Outline of the QCD Analysis 



In order to study the applicability of pQCD quantitatively, the (2+1) jet event rate for different 
bins of Q 2 (see table |l|) is chosen |3|, |l9| to determine ot s {ii 2 ) as a function of the renormalization 
scale^ /J%, which is taken to be Q 2 . The NLO QCD calculation can be used to calculate an 
expression for -R2+1 (Q 2 > Vcut) in terms of a power series in a s (Q 2 ) with coefficients a(Q 2 ,y cu t) 
and b(Q 2 ,y cut ): 

R2+i{Q\ycut) = a2+l{ ^ ycut \ = a(Q 2 ,y cut )a s (Q 2 ) + b(Q 2 ,y cut )a 2 (Q 2 ) + 0(a 3 s ). 

Taking the variation of a s with Q 2 to be negligible over each bin of Q 2 , to 0(a 2 ) this equation 
reduces to a quadratic form for a s with coefficients determined from the integrals of a(Q 2 ,y cu t) 
and b(Q 2 , y cu t) over each Q 2 bin. This quadratic expression can then be solved for a s using the 
measured jet event rates for different Q 2 . 

The cross section used for normalization, cr norm (Q 2 ), is taken as the sum of jet cross sections, 
where one or two current jets for (1+1) and (T2+1 f° r (2+1) events) are found in the 

acceptance region. Experimentally R2+1 is obtained as N2+1 / (N2+1+ N\ + i) , where Ni+\ denotes 
the number of events with % observed current jets. 

The theoretical 0(a 2 ) Monte Carlo integration programs, MEPJET fjj] and DISENT g], 
used for the prediction of jet cross sections are only available at the parton level and not as 
event generators. Both programs give the same results to within 1% precision. Consequently 
the relation between observed jets in the detector and the underlying partonic jet configuration 
has to be obtained in form of correction factors from a LO event generator which includes a 
simulation of higher-order effects and hadronization. 

The MEPJET Monte Carlo program is used to compare the NLO prediction with the mea- 
sured jet event rates and distributions of kinematic variables. It is necessary to apply to the NLO 
Monte Carlo the same phase space restrictions in terms of positron variables and jet quantities 
as is done in the experimental event selection. The forward direction is the most crucial region 
for the comparison. The emission of multi-gluon radiation from initial state partons entering the 
hard scattering process causes a higher reconstructed jet multiplicity. The influence of multi- 
gluon radiation increases with decreasing Q 2 because at low Q 2 very low values of the scaling 
variable xb = Q 2 /2qP are reached, which increases the permitted phase space (£ > xb) for the 
momentum fraction £ of the proton carried by the parton initiating the multi-gluon radiation. 
This initial state radiation is naturally not taken into account by a fixed order theory (as in 
MEPJET). As a consequence the jet event rate which is predicted by the NLO calculations is 
much smaller than the measured event rate without further cuts, particularly at low values of 
Q 2 . Therefore the event generator Monte Carlo is not only used to calculate the correction 
factors but also to determine cuts necessary to suppress the influence of multi-gluon radiation 
and hadronization effects. 

To define such a phase space region where the NLO QCD prediction can describe the partonic 
scattering underlying the data events, an event generator is necessary, which has a parton level 
comparable to the partonic configuration of the NLO QCD calculation and describes the data. 
Therefore the acceptance region for the current jets is chosen using the following criteria: 

• The dependence of the jet event rate on Q 2 predicted by the NLO calculation should be in 
good agreement with the prediction of the event generator at the parton level. For a rea- 
sonable range of the QCD parameter A^L as used in MEPJET^ (150 < A^L < 400 MeV) 

2 The factorization scale fi 2 - is also set to Q 2 . 

3 In MEPJET the effective number of flavors n/ contributing to the QCD parameter A. * (the calculation 
assumes zero quark masses) is fixed to 5; but in the text presented here always the corresponding AilL will be 
given. 
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also the absolute values of the jet event rate should be close to the predicted values from 
the event generator in each bin of Q 2 . In addition other reconstructed jet variable distri- 
butions at partonic level, such as jet angles, should be in acceptable agreement with the 
NLO calculation. 

• The measured (2+1) jet event rates and the other data distributions of the chosen jet 
variables should be described by the event generator model at the detector level. 



It is assumed that if these requirements are fulfilled the correction factors from measured jets 
to the NLO parton configuration can be provided by the event generator. 



As event generator the program LEPTO 6.5 [21] is used. LEPTO enables event generation by 



matching leading log parton showers to the LO matrix elements (MEPS). In the event generator 
a momentum scale Qq specifies the termination of the parton shower, that is the boundary 
between the perturbative and non-perturbative stages of jet development. For this analysis, this 
scale is taken to be 1 GeV as motivated in [22]. It has been demonstrated that this combination 



of matrix elements and parton showers gives a good description of the data in a wide range of the 
phase space specified by the cuts below (3|, 2C, 23, 24]. Nevertheless, there is some arbitrariness 
in the definition of the parton level. Investigations have shown that the determination of a s 



can depend on the definition of the parton level at the order of a few per cent [25, 26]. The 
interpretation of these investigations is, however, difficult. It can be argued that this sensitivity 
is due to effects of missing higher order terms which may well be included in the systematic 
errors due to the uncertainties of the available QCD models, in the renormalisation scale and in 
the choice of the used jet algorithm (see below). Accordingly, no additional source of systematic 
error is included for this arbitrariness. 

For the selected phase space region, in each Q 2 bin a correction factor to the experimental 
jet event rate is calculated using the LEPTO 6.5 generator from the ratio of the jet event rates 
obtained at the parton level after parton showering and at the detector level. In the same 
way also correction factors from the detector to the hadron level are determined. A possible 
model dependence introduced by this prescription is investigated by using the event generators 
HERWIG @ and ARIADNE |28j and its effect is included in the systematic error. 



The theoretical predictions are calculated using parton density parameterizations with the 
factorization scale /i 2 given by Q 2 . The parton densities are usually obtained via evolution from 
a low lying scale using a value a s (M|) as input parameter, which may differ from that extracted 
in the analysis presented here. It will be shown that the error introduced by this inconsistency 
is small in the kinematic region considered. 



6 Phase Space Selection 

The events selected according to section 3 are subject to a jet classification applying the mod- 
ified JADE jet algorithm as discussed in section ||. A resolution parameter of y cu t = 0.02, in 
combination with the cut W 2 > 5000 GeV 2 (section ||), implies an invariant mass squared of 
the two current jets which is greater than 100 GeV 2 . After the jet classification a cut in the 
maximum jet angle {9 max = 145°) is applied to ensure that the jets are measured in the LAr 
calorimeter. 

A cut at moderate values of Q 2 (Q 2 > 40 GeV 2 ) ensures that R2+1 is well described by 
LEPTO. In addition it avoids phase space regions where the variation of the QCD parameter A 
is more important in the parton density functions than in the hard-scattering cross sections 2S]. 



Furthermore, since this cut does not suppress parton showers sufficiently, additional requirements 
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are necessary. Parton showers produce jets with predominantly small angles 
the incident proton direction and small values of 



relative to 



mm 



Ej e t,i{^ 



cos 9 



jet,i) 



=1 ' 2 Ej=l EjetjO- ~ COsOjetj) 



The values of z p determined by the measured energies Ej et j and angles Ojetj (j = 1> 2) of the jets 
in the laboratory frame are in the range < z p < 0.5. The 9j e t and z p distributions are strongly 



correlated but are projections of different regions of the jet phase space in the 
Both are used in this analysis to restrict the jet acceptance region. 



"jet > Zp 



)-plane. 



In fig. |3)a the distribution of the smallest jet angle (9jet,mm) relative to the incident proton 
direction in a (2 + 1) jet event using the JADE algorithm is plotted for Q 2 > 40 GeV 2 and 
z p > 0.1. The curves are normalized to the number of events. It can be seen that above 
9jet,min = 10° the data are well reproduced by the MEPS model after detector simulation. In 
addition the MEPS model at the parton level and the MEPJET calculation (with MRSH as 
parton density parameterization and = 230 MeV) agree reasonably well with each other. 
Both the MEPS model and the NLO calculation predict however fewer (2 + 1) jet events with 
smallest 6j e t,min values than measured in the data. A 9j et cut suppresses parton showers and 
hadronization effects efficiently. The z p distribution is shown for Q 2 > 40 GeV 2 and 6j e t > 10° 
in fig. ||b. The final phase space selection fulfilling the criteria described in section || is therefore 
given by the requirements Q 2 > 40 GeV 2 , 6j e t > 10° and z p > 0.1, and is hereafter referred to 
as the central cut scenario. With these requirements the final event sample contains 1038 (2+1) 
jet events of which 689 have Q 2 > 100 GeV 2 . 



Due to effects of the finite resolution of energy and angles in the calculation of m 2 j , there are 
migrations between jet classes at the levels of the detector, hadrons and partons. The efficiency 
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Figure 3: a) The 9j e t t min distribution of the (2+1) jet events for Q 2 > 40 GeV 2 , z p > 0.1 
and 9 max = 145°; b) The z p distribution of the (2+1) jet events for Q 2 > 40 GeV 2 , 
10° < Ojet < 145 . The JADE algorithm is used to define the jets; in both figures the data 
are represented by the points with statistical error bars only. The solid line shows the prediction 
of the MEPS model at the detector level and the dashed line the prediction of the same model at 



the parton level. The NLO calculation (A- 



U)_ 

MS 



230 MeV) is given by the dotted line. All curves 



are normalized to the number of (2+1) jet events. 
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Figure 4: The a s determination at four values of Q and the evolution to the scale Q = M% 
using the modified JADE jet algorithm. The data points are placed at the average Q 2 in each 
bin. 

a) The measured R,2+i(Q 2 ) (points) and the jet event rates corrected to the parton level (tri- 
angles) together with the prediction of the MEPS Monte Carlo at the detector level including 
statistical error bands (solid line), parton level (dashed line) and the MEPJET NLO calcula- 
tions (dotted and dashed dotted lines) for two values of A^-. The event rates are determined 

for y cu t = 0.02, Q 2 > 40 GeV 2 , 6 ma x = 145° , 9j e t > 10° and z p > 0.1. The predictions are 
based on the MRSH parton density parameterization. 

b ) The a s values derived from the jet event rates as function of Q 2 with statistical errors from 
the data and the correction factors as inner error bars, and total errors as full error bars. The 
width of the Q 2 bins are indicated by the horizontal bars. The fit result (dashed line) and the 1 
s. d. errors ( dotted lines ) are based on the individual a s ( Q 2 ) values with their statistical errors 
only. The solid lines represent the obtained range for a s (M^) taking into account the systematic 
uncertainties, except the jet algorithm error. The open triangle at Q 2 = M| indicates the world 
average a s (M§) = 0.118 ± 0.003 JU/. 



and purity of the selected BEMC (LAr) sample are 27% (40%) and 57% (69%) with respect to 
the jet classification at the parton level and the above cuts p5| ]. 

The influence of radiative QED corrections on the (2+1) jet event rate is studied with the 
event generator DJANGO |3(| and is negligible. In this acceptance region the NLO corrections 
to the LO calculation are less than 30% in the investigated range of Q 2 and the Bjorken variable 
x B - 

Using the JADE algorithm the measured (2+1) jet event rate -R2+1 for Q 2 > 40 GeV 2 , 
9 jet > 10°, z p > 0.1 is shown in fig. ||a. The data above Q 2 = 40 GeV 2 are divided into four 
bins containing approximatively equal numbers of (2+1) jet events. The measurements compare 
reasonably with the MEPS Monte Carlo at the detector level. In fig. ||a it can also be seen 
that the jet event rate at the parton level in the MEPS model lies well between the MEPJET 
prediction for A^L = 150 MeV and A^L = 400 MeV. 
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Q 2 [GeV 2 ] 


R2+1 [lO" 3 ] 


R had [ 1Q -3] 


flTx [10- 3 ] 


f(det — > /iad) 


/(dei — > par) 


40-100 


25 ± l±l 


33 ±4l^ 


47 ± 6i | 


1.32 ±0.15 + 0.14 


1.88 + 0.21-0.03 


100-300 


46±2±^ 


60 + 4^ 


73 ± 5± 8 9 


1.30 ±0.04 +0.00 


1.59 + 0.05 - 0.12 


300-700 


57 ± 4± g 


72±5t^ 


91 ±61 J 


1.26 + 0.04- 0.02 


1.60 + 0.05 - 0.14 


700-4000 


69 ± 6+t 


82±8l^ 


99 ± 10±^ 


1.19 + 0.06-0.02 


1.43 + 0.07-0.09 



Table 1: The measured (2+1) jet event rates {R2+1) with statistical and systematic errors 
in different bins of Q 2 for the central cut scenario using the modified JADE algorithm and 
corrected to the hadron (R^'+i) an d parton level (R^j). In addition the correction factors from 
the detector to the hadron {f{det — > had)) and parton (f(det — > par)) level are included with 
errors due to Monte Carlo statistics and the model uncertainty estimated by ARIADNE. The 
statistical errors on the measured R2+1 are given by the data only, whereas the errors on the 
corrected (2+1) jet event rates (R^^, R V 2+i) include also the statistical errors from the Monte 
Carlo. The systematic errors on R2+1 correspond to the uncertainty on the hadronic energy 
scale. The systematic errors on ^ take into account the hadronic energy scale and model 
uncertainty. For R^+i also the effects of parton showers, hadronization parameters and QCD 
Monte Carlo models were taken into account similar to the error estimation on a s as described 
later in the text. 



7 Determination of the Strong Coupling Constant 

The correction factors f(det — > par), necessary to obtain from the measured jet event rate R2+1 
the rate at the parton level (-Rf+l)' are gi ven i n table |] for the different Q 2 bins shown in fig. |4]a. 
They are between 1.4 and 1.9. Also the factors f{det — > had) for correcting the measured event 
rates to the hadron level are included in table [l|. It can be seen that hadronization and detector 
effects lead to lower observable (2+1) jet event rates than at the parton level due to resolution 
effects. The statistical errors on the corrected event rates in table [l] include the statistical 
error from the data as well as from the calculation of the correction factors using the Monte 
Carlo generator. For all calculations the parton density distributions were parameterized using 
MRSH |2|. 

The experimental event rates corrected to the parton level i^f+i are a ^ so plotted in fig. |||a for 
the standard set of cut parameters (see above) using the JADE algorithm. They are converted 
into four values of a s (Q 2 ). The QCD parameter A^L in the 2-loop solution of the renormaliza- 
tion group equation is fitted [|^] to these extracted values of a s (Q 2 ) taking into account their 
statistical errors only. The results, together with the ±1 s.d. parameter error of the fit, are 
shown in figure |]b. The observed Q 2 dependence is compatible with the expected theoretical 
evolution. The numerical result at Q 2 = Mf is a s (M§) = 0.117 ± 0.003 (stat). 

To obtain the systematic uncertainty on a s (M§) the fit is repeated for various other scenarios 
changing cuts and parameters relevant for the jet analysis. This procedure automatically takes 
into account the correlation between the four bins in Q 2 . All systematic errors discussed in 
detail below are summarized in table 

In order to determine the systematic experimental error the following kinematic cuts are 
investigated: 

• The cut in z p is varied between 0.05 and 0.15 for 8j e t > 5°, 10°, 15° leaving the other 
parameters at their central values. In addition 9 max is varied between 140° and 150° 
leaving the other parameters at their central values. The resulting error on a s (M z ) as 
determined from the maximum spread of all fit results is +0.004 and —0.004. 
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• The hadronic energy scale of the LAr calorimeter in the data reconstruction is varied by 
±4%. The resulting error on a s (Mf ) is +0.005 and -0.003. 



The systematic experimental error (+0.006 and —0.005) on a s (M|) is taken as the quadratic 
sum of the two errors given above. Effects due to the uncertainty of the energy of the scattered 
positron and the variation of kinematic cuts defined in section 3 are negligible. 

A further source of systematic uncertainty stems from the dependence of the correction 
factors on the Monte Carlo models used to calculate them. These uncertainties are analyzed 
employing the following procedure: 



The LEPTO Monte Carlo uses the JETSET [34] fragmentation model. Several sets of 



parameters describing the hadronization and fragmentation as fitted by the LEP experi- 



ments [35] are used to study the stability of the correction factors from the hadron to 
the parton level. The widths of the Gaussian distributions for the transverse momenta 
of partons in the proton (kx (PARL(14) ) and a p t (PARL(3) )) are varied between 350 — 
700 MeV and 440 — 700 MeV, respectively; the fragmentation parameter a of the LUND 
fragmentation functions is varied between 0.1 and 1; the width a (PARJ(21)) in the 
gaussian p x and p y transverse momentum distributions for primary hadrons is changed 
from the default value (360 MeV) to 250 and 450 MeV. In addition the momentum scale 
which determines the termination of the parton shower (Qo-P arame ters PYPA(22) and 
PARJ(82) in LEPTO for the initial and the final state partons respectively) are varied 
between 1 and 2 GeV and between 1 and 4 GeV, respectively. Finally the lowest value 
of the invariant mass squared s of the two generated partons, for which the LO ME 
calculations are carried out, is changed from 4 GeV 2 to 25 GeV 2 . The resulting fitted 
a s (M|) values are in the range 0.115 to 0.121. 

Instead of LEPTO, the HERWIG @ Monte Carlo model with its different treatment 
of parton showering and hadronization is applied to calculate the (2+1) jet event rate 
on parton and hadron level. HERWIG is able to describe the shape of the measured 
i?2+i(Q 2 ) distribution. Using the same corrections for detector effects as determined with 
the LEPTO generator, the fit leads to a reduction of a s (M|) by 0.007. 



Source 


Variations 


Aa s (M 2 ) 


Statistics 




±0.003 


9 j e t and z p cut 


variation in 6j e t and z p , see text 


+ 0.004 

-0.004 


Hadronic energy scale 


±4% 


+ 0.005 
-0.003 


Model dependence 


LEPTO 6.5 


+ 0.004 

-0.002 


HERWIG 5.9 


-0.007 


ARIADNE 4.08 


-0.005 


Parton density functions, 
used in the PDF's 

MS 


MRSAp-201, GRV92, CTEQ-4A1 
0.15,0.40 GeV 


+ 0.003 
- 0.002 


Renormalization and factorization scale 


/i? = l/4,4Q 2 i=f,r 


+ 0.003 
-0.007 


Vcut 


0.015, 0.025 


+ 0.002 
-0.003 


jet algorithms 


E0 
P 


+ 0.000 
+ 0.006 



Table 2: The determination of various systematic uncertainties. The central cut scenario and 
the MRSH parton density parameterization are used as reference. 
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• The event generator ARIADNE 4.08 |2f| is also used to study the model dependence. 
This Monte Carlo model treats the gluon emission in the framework of the colour dipole 
model [36] as an alternative to the parton showers in LEPTO. In the selected phase space 
ARIADNE describes the data well with exception of the i?2+i(Q 2 ) distribution. The fitted 
a s (M|) value using fully simulated ARIADNE events deviates by —0.005 from the quoted 
result obtained with LEPTO. 



The dependence of a s (M§) on the choice of the Monte Carlo model is estimated by the quadratic 
sum of all these uncertainties to be +0.004 and —0.009. 

Finally the theoretical uncertainty is investigated by changing various input conditions for 
the NLO calculations: 



The MRSH parton density function (PDF) used in MEPJET is replaced by the MRSAp- 
201 @, GRV92 |H, and the CTEQ-4A1 || parameterizations. They all use a similar 



A^L value of about 230 MeV. The largest change in the fitted value of a s (M§) is ±0.001. 



As a particular point of interest 1 4C , |28|] the correlation of the measured value of a s in this 



analysis with the value used while fitting the PDFs is investigated. This is done using the 



GRV92 PDFs which are available for various A^L values El in the range 0.15-0.40 GeV. 

MS 1 — V ° 

Although a s (M§) used in the PDFs changes by up to 0.017, the fitted values only differ by 
+0.003 and —0.002. This large reduction is due to the fact that the value of a s influences 
the scaling violation, whereas i?2+i(Q 2 ) is directly proportional to the strong coupling 
constant. In table [2] both errors were combined to one. 

The factorization scale fi 2 and the renormalization scale fi 2 are varied independently be- 
tween Q 2 /A and AQ 2 . The effect on the measured value of a s (M|) is +0.003 and -0.007. 

The dependence on the chosen value of y cu t was investigated by varying y cut between 0.015 
and 0.025 in the NLO calculation, LEPTO simulation and for the data. Above 0.025 the 
number of (2+1) jet events rapidly decreases but the results are still consistent within the 
larger statistical error. Below 0.015 the regime of very small jet masses opens up, making 
the fit procedure unstable. The resulting error on a s (M z ) is +0.002 and —0.003. 



The total theoretical uncertainty (+0.005 and —0.008) on a s (M|) is obtained summing these 
three errors in quadrature. 

The result of the analysis using the JADE algorithm, where the total systematic uncertainty 
is taken as the quadratic sum of all uncertainties discussed above, is: 

a a (M§) = 0.117 ± 0.003 (stat) t (sys) . 

The solid lines in figure |4|b indicate the full range of a s (Q 2 ) obtained in this experiment using the 
JADE jet algorithm adding statistical and systematic errors in quadrature at the scale Q 2 = M§. 
They correspond to a s (M§) = 0.104 and 0.126. 

To obtain the systematic uncertainties on the four measured values of a s (Q 2 ) using the JADE 
algorithm the following method is applied. Each of the used scenarios, except the variation of 
the renormalization scale, results in measured values of a s (Q 2 ) for the four bins in Q 2 . From the 
values of the different cut scenarios the systematic experimental and theoretical uncertainties, 
as well as the model dependence in each bin, are determined in the same way as discussed above 
for the uncertainties on a s (M§). The total systematic uncertainty using the JADE algorithm is 
given by the quadratic sum of these three error contributions. The outer error bars of the a s (Q 2 ) 
points in figure ||b represent the total uncertainties, which are calculated as the quadratic sum 
of statistical errors and total systematic uncertainties (see table ||) . 
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JADE 


E0 


P 


Q 2 [GeV 2 } 


a s (Q 2 ) 


stat 


sys 


tot 


a s (Q 2 ) 


a s (Q 2 ) 


54 


0.184 


±0.017 


+ 0.035 
- 0.038 


+ 0.039 
- 0.042 


0.189 


0.206 


209 


0.166 


±0.009 


+ 0.028 

-o!o2i 


+ 0.029 
-o!o23 


0.166 


0.177 


440 


0.151 


±0.010 


+ 0.018 
- 0.022 


+ 0.021 
- 0.024 


0.151 


0.167 


1272 


0.130 


±0.012 


+ 0.012 
- 0.020 


+ 0.018 
- 0.024 


0.127 


0.128 


Ml 


0.117 


±0.003 


+ 0.009 
-0.013 


+ 0.009 
-0.013 


0.117 


0.123 



Table 3: The values of a s (Q 2 ) together with their statistical, systematic and total uncertainties 
for the modified JADE algorithm at different values of Q 2 , which are taken to be the average 
in each considered Q 2 range (see table 1). In addition the extracted a s (Q 2 ) values using the 
modified E0 and P algorithms are included. 



To study the dependence of the determined a s (M z ) value on different jet algorithms the 
analysis is redone for the central cut scenario using the E0 and P variants of the JADE jet 
algorithm. The obtained a s values in the different bins of Q 2 are given in table ^. The extracted 
value of a s (M z ) based on the E0 algorithm, a s (M z ) = 0.117, is, as expected from a theoretical 
point of view, very close to the value determined using the JADE algorithm, whereas for the P 
algorithm the fitted value is increased by ±0.006. The final result of this analysis is therefore 

a s {M 2 z ) = 0.117 ± 0.003 (stat) ± (sys) ± 0.006 [jet algorithm) , 

where the last error reflects a possible sensitivity to the choice of the jet algorithm used. 

Another approach taken by the HI Collaboration is to perform a quantitative pQCD analysis 



based on jet event rates differential [26| in the JADE jet resolution parameter. This method using 
a similar phase space region leads to a s (M§) = 0.118 ± 0.002 (stat) ILq'qqI ( s U s ) io'oofl (theory), 
fully consistent in value and magnitude with the result of the present analysis. 

The combined result from studies of the hadronic final state in the experiments at LEP and 
SLC |3lJ is a s (M%) = 0.122 ± 0.007. Here the error is totally dominated by theoretical and 
model dependencies associated with the same error sources as in the analysis presented in this 
paper. 



8 Summary 

Jet production in neutral current deep inelastic ep scattering at HERA is studied using the 
modified JADE jet algorithm. The strong coupling constant a s is determined over a wide range 
of Q 2 by evaluating jet event rates corrected to the parton level assuming the hadron-parton 
relation as given by the event generator LEPTO. The Q 2 dependence of the determined a s 
values is compatible with the theoretical prediction of the renormalization group equation. A fit 
of the QCD parameter A^L leads to a s (M 2 z ) = 0.117 ± 0.003 (stat) ± Sloxs ( s V s ) ■ The systematic 
error is dominated by the uncertainty in the hadronic energy scale, the renormalisation scale 
dependence, and the dependence on the Monte Carlo model. Adding statistical and systematic 
errors in quadrature yields a s (M z ) = 0.1171 o!o?3 consistent with the world average a s (M z ) = 



0.118 ± 0.003 [31]. A dependence of the extracted a s (M z ) value on variations of the JADE 



algorithm is found. The largest change observed is ±0.006, leading to the final result: 
a s (M%) = 0.117 ± 0.003 (stat) ± (sys) ± 0.006 (jet algorithm) . 
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